Quantum annealing and the Schrodinger-Langevin-Kostin equation 
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We show, in the context of quantum combinatorial optimization, or quantum annealing, how the 
nonlinear Schrodinger-Langevin-Kostin equation can dynamically drive the system toward its ground 
state. We illustrate, moreover, how a frictional force of Kostin type can prevent the appearance 
of genuinely quantum problems such as Bloch oscillations and Anderson localization which would 
hinder an exhaustive search. 
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I. INTRODUCTION 

The aim of combinatorial optimization is to find good 
approximations to the solution(s) of minimization prob- 
lems. Many of the most famous algorithms currently used 
in this field [T were inspired by analogies with physical 
systems. Among them the most celebrated is Thermal 
simulated annealing |2] proposed in 1983 by Kirkpatrick 
et al.: the space of all admissible solutions is endowed 
with a potential profile dependent on the cost function 
associated to the optimization problem. The exploration 
of this space is represented by a temperature dependent 
random walk. An opportunely scheduled temperature 
lowering (annealing) stabilizes then the walk around a, 
hopefully global, minimum of the potential profile. 
The Quantum annealing approach to combinatorial op- 
timization [2111], instead, was originally suggested by the 
behaviour the stochastic process q,j{t) associated to 
the ground state of a Hamiltonian of the form: 

where the potential function V encodes the cost function 
to be minimized. The behavior of qi, is characterized by 
long sojourns around the stable configurations, i.e. min- 
ima of V{x), interrupted by rare large fluctuations which 
carry q^, from one minimum to another: q^, in thus al- 
lowed to "tunnel" away from local minima to the global 
minimum of ^(x). The deep analysis of the semiclassical 
limit performed in [7 shows, indeed, that as z/ 0+ "the 
process will behave much like a Markov chain whose state 
space is discrete and given by the stable configurations" . 
However, the ground state of H^, is seldom exactly known 
and approximations are required. One of the earliest pro- 
posals in this direction, advanced in [4 and applied in 
[3 , was to construct an unnormalized approximation of 
(/>j/(a;) by acting on a suitably chosen initial condition 
4>triai{x) with the Hamiltonian semigroup exp(— tiJ^), 
namely by solving, with the initial condition (ptriai, the 
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imaginary time Schrodinger equation. Similar ideas ap- 
pear in the chemical physics literature [8] and, with 
more specific reference to the optimization problems con- 
sidered here, in [51 UHl lH]- Yet, the inability to au- 
tonomously construct the ground state process, without 
recourse to the unphysical step of imaginary time evo- 
lution, substantially detracts from what is otherwise a 
physical route to optimization by dynamical evolution 
toward the ground state. 

In this note, encouraged by the progress in quantum an- 
nealing in the last twenty years, as reviewed for instance 
in [121 [131 (Hj, by its close relationship with adiabatic 
quantum computation [15] and by proposals of its hard- 
ware implementation |16j . we try to eliminate the above 
unphysical step: we try to implement, instead of imagi- 
nary time evolution, the idea of reaching the ground state 
with the help of viscous friction [TT] [HI [19] . 

We first introduce the nonlinear Schrodinger-Langevin- 
Kostin (SLK) equation and illustrate, by means of ex- 
amples on two toy models, how a frictional force acts in 
the continuous case. Then we turn to quantum combina- 
torial optimization and show how dissipation can, in the 
discrete case, balance genuinely quantum effects, such as 
Bloch oscillations and Anderson localization, which can 
hinder the search of optimal solutions. 

II. CONTINUOUS CASE 

The SLK equation is the analogue of the Heisenberg- 
Langevin equation and represents a quantum analogue 
of classical motion with frictional force proportional to 
velocity [TT] |TH]; it can be seen as a rough analogue of 
the classical Drude-Lorentz model of Ohmic friction, i.e. 
an approximate description of the motion of a quantum 
particle through matter with inelastic scattering. 
A solution ip{t,x) = \J p(t, x)e'"^(*''^^ of the equation: 

.dibit, x) v"^ d'^ib{t,x) ,w ^ 

+ 13 S{t,x)il^{t,x) (2) 

satisfies the inequality ^( V-C^) l^^i^l V'(i) ) < for /3 > 0, 
H,j being the Hamiltonian ([T]). What we will show be- 
low is how the norm preserving dissipative evolution de- 
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FIG. 1: (Color online) Frame (a): The probability density 
a;)p at the initial time t — and final t = tmax ~ 50. 
The potential V{x) is drawn in arbitrary scale for expository 
purposes. Frame (b): The expected value { tf}{t) |_ff^| ^p{t) ) of 
the Hamiltonian operator (i/ = 1, f3 — 0.5); in the inset: the 
overlap |{ i/>(t) | )|^ between the ground state oi Hi, and the 
state of the system at time t, for < t < tmax = 50. 



scribed by (pi) can dynamically drive a suitable initial 
condition -0(0,2:) toward the ground state (j>,^{x) of H^^. 
Toy model 1: Require 



(/)^{x) = c+ exp 



{x — aY 



c_ exp 



{x + a) 

4(72 



Co exp 



(the parameters c±^o being chosen so that 4>v{x) > 0) 
to be the ground state of H^, and to belong to the 
eigenvalue E,^ = 0. The above two requirements de- 
termine the potential V{x) = (x) ^'^'^(^)- fig- 
ure [l] follows the evolution ip{t, x) of the initial condi- 

under ( 2 1 for a time interval 



tion "0(0, x) 



It shows that, as {ip{t) \Hj\y;{t)) decreases with 
time, the "vacuum overlap" \{ip(t) \ (f>^ )\ approaches the 
value 1. While the state ip{t) approaches the ground 
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FIG. 2: (Color online) Frame (a): The probability density 
\tl>{t, x)\^ at the initial time t = and final t = tmax = 50. The 
potential V{x), corresponding to the choice a+ — 2.25, a_ = 
1.75, Vo ~ 1, S = 0.1 of the parameters of ([3|, is drawn 
in arbitrary scale for expository purposes. Frame (b): The 
expected value ( tp{t) \H,^\ ^{^{t) ) of the Hamiltonian operator 
{v = 1, 13 = 0.3); in the inset: the potential profile W(t,x) at 
the initial and final times: the potential profile at t = tmax = 
50 completely overalps V{x). 



state, some of the probability mass "tunnels" from the 
leftmost (local) to the rightmost (global) minimum. 
We point out that this class of examples, where both the 
ground state wave function (/)y(a;) and the ground state 
energy E^^ are known, allows also for the calibration of 
the numerical method. In our case we have used the 
built-in NDSolve resource of Wolfram Mathematica 6. 

Toy model 2: For the sake of comparison with classi- 
cal literature on quantum annealing [12l [20] , we consider 
here a double-well potential of the form 

Vo^ 4^ + <5x, for X < 0. 



As shown in figure 2(a)| for the parameters used there, 
the local minimum of the potential V{x) is wider than 
the global one. We refer the reader to section 2.2 of [12] 
and to [20] for a discussion of the meaning of the param- 
eters and for the presentation of numerical experiments 
comparable with ours. Here, we use this well known toy 



model as an example in which the ground state is un- 
known and the dissipative dynamics of SLK type pro- 
vides a method to find it. 



While the initial condition evolves (figure [2(a) ), the mean 
value of decreases as in figure [2(b)] That ij{tmax,x) 
is a good approximation of the ground state is shown 
by comparing, in the inset of figure [2(b)[ V{x) with the 

potential W{t,x) = 2^ 2:) + ( V(i) l-^.l ^ ) 

evaluated at t = tmax, of which ipi'^maxi x) is the ground 
state belonging to the eigenvalue 0. Comparison of the 
two curves in the inset of figure [2(b)[ is, furthermore, 
suggestive of a real-time version of Piela's method of de- 
formation of the potential energy hypersurface [9\. 
In the continuous case, then, the SLK Hamiltonian en- 
riches quantum annealing of what can be seen as a proba- 
bility percolation: while the state of the system converges 
toward the ground state, some probability mass tunnels 
from one minimum to another and, instead of tunneling 
back, as would happen in a reversible dynamics, remains 
there. This autonomous stabilization would represent an 
alternative to adiabatic quantum computation |15j . 



III. DISCRETE CASE 

Most of the effort in Ref. [4J went into making the 
intuition developed so far available in a context of com- 
binatorial optimization. In such a context the domain 
of the function V to be minimized is a finite set Q and 
the search of the minimum of V is modeled on a graph 
{Q,E), where the edges e G E describe the moves al- 
lowed in the search. For instance, in [3], Q was taken 
to be the Boolean hypercube Qn = {— Ij l}", for some 
positive integer n, and an edge was placed between any 
two points in separated by a unit Hamming distance. 
In this note, we consider the much simpler instance 
in which Q is, for some positive integer s, the fi- 
nite set As = {l,2,...,s} equipped with the set 
of edges E ^ {{i,j} : £ A, x A, A \i - j\ = 1}. Ac- 
cording to the general approach outhned in [3], this 
amounts to a search of the minimum of the function V 
defined on A^ by means of an interacting continuous time 
quantum walk |21j on A^ governed by a Hamiltonian of 
the form 



h=- 



|j + i)(j| + b-)(.7 + i|+ 



nj-)ij>(j 



(4) 



We wish to show here that the quantum search outlined 
above can suffer from two, typically quantum, problems, 
namely Bloch oscillations [S^ and Anderson localization 
[53] and that a certain amount of "viscous" friction can 
provide some relief to both these problems. 
On a finite box As, we consider the evolution of an ini- 
tial condition of the form Ck{x) = X)^= 
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FIG. 3: (Color online) s = 100, e = 17, fc = (e - l)/2, go = 
2/s, < t < 4s. Frame (a): free evolution (g = 0, (3 = 0); 
(b): g — 3go, /3 = 0: Bloch oscillations prevent the wave 
packet from exploring most of the solution space; (c): g = 
3go, 13 = 4go: Bloch oscillations no longer appear, thus al- 
lowing a complete exploration of the solution space and con- 
vergence to the ground state; (d) the probabilies, as a function 
of time, of reaching the 5 = 2e rightmost sites of As in the 
free case (thin dashed line) and in the forced and damped 
case g = 3go, fH = 4go (solid thick line). 



We refer the reader to section 5 of [2T[ for a motiva- 
tion of this choice: suffice here to say that it describes 
a spatially well located wave packet that in the absence 
of any potential moves back and forth, with speed close 
to 1, inside the box As, as in figure 3(a)[ The effect on 
this ballistic evolution of a linear potential V{x) = —gx 
is shown in figure [3(b)[ For g = 0(l/s), the peculiar 
energy-momentum relation E(p) = 1 — cosp, holding on 
a discrete lattice, determines Bloch oscillations that pre- 
vent the wave packet from approaching the point x = s 
at which the minimum of the cost function V{x) is lo- 
cated. Figure [3(b)| is therefore a reminder of the fact 
that a greedy quantum optimization driven by the cost 
function itself acting as a potential can be hindered by 
the fact that, on a lattice, increasing momentum p can 
mean decreasing velocity v{p) — sin p. 

We propose here to introduce a certain amount of vis- 
cous friction in the discrete Schrodinger equation, as a 
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"Kostin potential" K{t, x) (3 S{t, x) = (3 Arg {ipit, x)): 
.dip{t, x) 



dt 



{hiP){t,x)+K{t,x)ij{t,x)= (5) 
-l-{^{t,x + l) + ^t,x-l)) + 



+ V{x)'ip{t, x) + K{t, x)ilj{t, x). 

The idea is that friction can prevent the momentum p 
from crossing the first Brillouin zone and thus can pre- 
vent the velocity sinp from being inverted before the 
wave packet reaches the boundary of As . This unwanted 
inversion is illustrated in figure |3(b)[ the effect on it of a 



suitable Kostin potential is shown in figure 3(c) 
As it is easy to check that, for ipjTj evolving 
according to (|5]), it is 



dt 



; m \h\ m ) = - E v/R^+^M^- 



(6) 



• {K{t, a; + 1) - K{t, x)) sin {S{t, a; + 1) - S{t, x)) , 

the actual form of K{t, x) that we adopt in or- 
der to achieve decrease of ( ^{t) \h\ 'ip{t) ) is 
K{t, x) = PY.l=2 sin iS{t, y) ~ S{t, y-1)), with /? > 0. 

Figure 4(a)| shows, instead, the effect, in the form of 
Anderson localization, of a random Gaussian potential 
of mean and variance a^, acting independently on 
each site of A^. The order of magnitude ctq = 
of the noise parameter a is suggested by a scaling 
argument [21] . 

Whereas the fact that friction can wipe out Bloch 
oscillations is well known [H], the less well known fact 
that we show here is that the pseudo-ballistic motion 
shown in figure 3(c)| is much more stable than the truly 
inertial motion represented in figure 3(a) with respect to 
the onset of Anderson localization. 



IV. CONCLUSIONS AND OUTLOOK 

As a final remark, we observe that the same framework 
developed in this note for the combinatorial optimization 
metaphor can be used, with minor changes, to describe 
an excitation travelling along a spin chain or a light pulse 
propagating through a waveguide lattice [25]. We conjec- 
ture, therefore, that SLK dynamics can be exploited also 



in those fields. For example, we can, maybe, increase 
the fidelity of state transmission, in presence of imper- 
fections, along a spin chain, by applying a "tension" at 
both ends of it [55] (see figure [4(d)| . The sole convergence 
toward the ground state could, instead, find applications 
in all-optical switching of light in waveguide arrays |27j : 
the injected light pulse can be steered toward a given po- 
sition by a suitable tuning of the thermal gradient which 
determines the potential profile of the lattice. Future 
work should be devoted to further investigation of this 
open research problems. 
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FIG. 4: (Color online) s = 100, e = 17, fc = (e - l)/2, go = 
2/s, o-o = (10/s)^''^ , < t < 20s. Frame (a): g = 0, /3 = 
0, a — 2(70 : the wave packet gets confined in the first half of 
As; (b): for <; = 0, (3 = 0, a = 2ao the probability of ever 
reaching the 5 — 2e rightmost sites is negligible (Anderson 
localization); (c): g — 3go, /3 — 'lgo,a = 2ao: viscous friction 
allows the particle to drift to the right, by successive sojourns 
(the vertical strips) around successive minima of the Anderson 
potential; the ensuing slow transfer of the probability mass to 
the right of As is shown in frame (d) . 
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